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We present improved measurements of the branching fractions for the decays -S" — > Df7v~ 
and B° -^ D^ K~ using a data sample of 657 x 10^ BB events collected at the T{AS) reso- 
nance with the Belle detector at the KEKB asymmetric-energy e'^e~ collider. The results are 
B {B° -^ D+-K-) = (1.99 ± 0.26 ± 0.18) x lO"'^ and B {B° -^ D+K') = (1.91 ± 0.24 ± 0.17) x 10-^ 
where the uncertainties are statistical and systematic, respectively. Based on these results, we de- 
termine the ratio between amplitudes of the doubly Cabibbo suppressed decay i?" — >■ D'^'k~ and the 



Cabibbo favored decay B" -^ D tt"*", Rdtt 
the last term denotes the theory error. 



= [1.71 ± O.ll(stat) ± 0.09(syst) ± 0.02(theo)] %, where 



PACS numbers: 13.25.Hw, 12.15.Hh. ll.30.Er 



In the standard model (SM), CP violation occurs due 
to a single irreducible phase appearing in the quark-flavor 
mixing matrix, called the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix p^, vifhich relates the weak interaction 
eigenstates to those of mass. Unitarity of the CKM ma- 
trix yields relationships between its elements that can be 
depicted as triangles in the complex plane [2']. B me- 
son decays offer a variety of ways to measure the angles 
and sides of the unitarity triangle (UT), formed from ele- 
ments in the first and third columns of the CKM matrix, 
and, hence to verify the CP violation mechanism of the 
SM. 

Of particular interest is the decay 5° -^ D^n^, which 
is dominated by the tree level b ^ u transition shown in 
Fig.[lja). Assuming SU(3) flavor symmetry, one can use 
this decay channel to determine the ratio between am- 
plitudes of the doubly Cabibbo suppressed decay B° — >• 
_D+7r~ and the Cabibbo favored decay iS" — >■ D~tt^ [3( 



Rdt! = tan 9c - — 




B{B" ^ D-7T+)' 



(1) 



where 9c is the Cabibbo angle, foifoj is the D{Ds) 
meson decay constant, and B{B'^ -> D+7r~) and B[B'^ — > 
i?~7r+) are the branching fractions of B^ — > Dfir^ and 
S" — > D^7r+. In Figs. [TJb) andlljc) we show the dom- 
inant Fcynman diagrams for the decays S" — > D^n^. 



The ratio Rdtt is an important input for the determina- 
tion of the UT angle 03, since the measurement of time- 
dependent CP violation in i?" — > D^tt^ [J| determines 
only the quantity Rdtt sm{2(j)i + 03)/(l + Rdtt)' where 
01 is the most precisely measured angle of the UT [^, [g] . 
Furthermore, it has also been suggested [7!| that the CKM 
matrix element \Vub\ (related to one of the sides of the 
UT) can be extracted from the measured branching frac- 
tion of B° -^ D+TT-. 

The decay i?° -^ DfK~ occurs via the internal W- 
exchange diagram [see Fig. [Hd)]. Potential contribu- 
tions arising from rescattering effects_l8[ could enhance its 
branching fraction. Recent studies tS], however, find the 
rescattering contribution to be negligible. Furthermore, 
the calculation of Ron in Ref. [3| assumes that the size 
of the VF-exchange amplitude in B^ — >■ D^tt^ [Figs.[IJe) 
and Wis)] is small compared to the corresponding tree 
amplitude. One can verify this hypothesis with an accu- 
rate measurement oi B{B'^ -^ DfK~). In the absence of 
rescattering the exchange diagram is the sole contributor 
to B^ — > DfK^, and hence it provides a measure of the 
W^-exchange contribution in B° — ^ D^n^. 

The decay channels B° -^ Dt^^^ and B'^ — >• DfK^ 
have been previously studied by the Belle [10| and 
BaBar [11| collaborations. In this paper, we present im- 
proved measurements of the branching fractions B{B'^ — > 




FIG. 1: Feynman diagrams for (a) B° -^ D^tt"; (b) the 
Cabibbo favored decay 5° -^ D~Tr'^; (c) the doubly Cabibbo 
suppressed decay B" — >■ D'^-k~; and the color suppressed W- 
exchange contribution to (d) B" — >■ Dt K~ , (e) S" — >■ _D~7r+, 
and (f) B° ->D+7r-. 



D+7r~) and ,B(i?° — ?► DfK~) based on a data sample of 
657 X 10^ BB decays, which is close to 8 times the size of 
the one used in our earlier result [lOJ . The data were col- 
lected with the Belle detector at the KEKB asymmetric- 
energy e+e" collider [l2|. A detailed description of the 
Belle detector can be found elsewhere [l3|. 

We select B° -^ D+n- and B° -s> D+K' decay can- 



didates [IJ] from events that have four or more charged 
tracks. Each track is required to be well measured in a 
tracking system that consists of a silicon vertex detector 
and a central drift chamber (both operating in a 1.5 T 
magnetic field), and to originate from the interaction 
point (IP). Track candidates must have a minimum trans- 
verse momentum of 100 MeV/c, and a distance of closest 
approach with respect to the IP less than 0.2 cm in the 
r-(j) plane, which is perpendicular to the z axis, and less 
than 4.0 cm along the z axis, where the z axis is defined 
by the direction opposite to the e"*" beam. Charged pions 
and kaons are identified by combining particle identifica- 
tion (PID) information obtained with various subdetec- 
tors: ionization energy loss from the drift chamber, time- 
of-flight information from an array of scintillation coun- 
ters, and the number of photoelectrons from an aerogel 
Cherenkov counter system. We distinguish kaons from 
pions using a likelihood ratio, TZk/v — ^^k/O-'K + Ltt), 
where Lk (L^) is the likelihood value for the kaon (pion) 
hypothesis. We require TZx/tv to be greater than 0.6 for 
kaon candidates, while tracks failing this requirement are 
classified as pions. The efficiency for kaon (pion) iden- 
tification ranges between 84% to 98% (92% to 94%) de- 
pending on the track momentum with a pion (kaon) fake 
rate of about 8% (16%). 



We reconstruct Df mesons in three decay modes: 
D+ -^ cj)TT+, K*°K+, and K°K+. The (j) {K*°) mesons 
are formed from K^K^ (K^n^) pairs having invariant 
masses that lie within 14 MeV/c^ (75 MeV/c^) of the nom- 
inal (J) (K*^) mass [l5|. Note that for kaons originating 
from a (j) decay we relax the TZk/-w requirement to 0.1 
due to the small background contribution. To reduce 
combinatorial background, we require jcos^if | > 0.3 for 
the D+ -> (/)7r+ {D+ -^ K*°K+) mode, where Or is 
the angle between decay products of the (j) {K*'^) and 
the flight direction of the Z?+ meson in the rest frame 
of the (j) {K*^). We reconstruct K^ mesons through 
the channel K'^ — ^ 7r"'"7r~, where we require the invari- 
ant mass of two oppositely charged tracks (with the 
pion mass hypothesis assumed) to be within lOMeV/c^ 
of the nominal K'^ mass [15[. The K'^ candidates 
must also satisfy momentum-dependent selection crite- 
ria based on their vertex topology and flight length in 
the r-(f> plane J16J . We select D^ mesons in a wide mass 
window (1.92GeV/c2 < Af^+ < 2.02GeV/c2), common 
to the three decay modes, for further studies. Finally we 
combine each D^ candidate with an oppositely charged 
pion or kaon to form a neutral B meson. 

For the reconstruction of B candidates we utilize two 
kinematic variables: the center-of-mass (CM) energy dif- 
ference, Ai? = Eb — fi'bcam, and the beam-constrained 
mass, Mbc = V^^bcam ~-Pb' where E'bcam is the beam 
energy, Eb and pB are the energy and momentum of 
the B candidate measured in the CM frame, respec- 
tively (c = 1 is assumed). The Afbc distribution for 
signal events peaks near the B mass, while the IS.E dis- 
tribution peaks at zero. We retain B candidates with 
5.27GeV/c2 < Mbc < 5.30GeV/c2 and -0.1 GcV < 
AE < 0.2 GeV. An asymmetric AE requirement is im- 
posed to suppress background contributions from B de- 
cays, such as i?" -> -D*+7r^ and B° -^ D^^K^ , at neg- 
ative /S.E values. 

About 5% of the selected B° -^ -D+tt" and B^ -^ 
DfK^ events contain multiple B candidates. In such 
cases we choose the one with the Mbc value closest to the 
nominal i?*^ mass 151. In order to determine the back- 



ground reduction criteria (described below), signal and 
background yields are estimated in the signal region with 
the requirements that \AE\ < BOMeV and that M^+ be 
within IBMeV/c^, 15MeV/c^ and 17MeV/c2 of the nom- 
inal Df mass for 07r+, K*'^K~^, and K^K"^, respectively. 
These requirements roughly correspond to a ±3cr window 
in resolution. 

Continuum e+e^ ^ <M {q = u, d, s, and c quarks) 
events are the dominant background. To discriminate the 
jet-like continuum background from signal we use mod- 
ified Fox- Wolfram moments [13 that are combined into 
a Fisher discriminant. We further combine the Fisher 
output with the cosine of the angle between the B flight 
direction in the CM frame and the z axis, to form a like- 
lihood ratio TZ = hgig/ (Lsig+hqq). Here, Lgig and hqq are 



the likelihood functions for signal and continuum events 
obtained with Monte Carlo (MC) simulations 18|. We 
impose separate requirements on TZ for the three decay 
modes in both B° -J> -D+tt" and B° -J> D+R-. These re- 
quirements are obtained by maximizing a figure of merit, 
S/y/S + B, where S and B are the number of signal and 
qq events expected in the signal region, calculated using 
MC simulated events. The requirements on TZ remove 
92% (78%) of continuum background while retaining 75% 
(86%) of signal events for 5° -^ D+tt" (B° -^ DfR-). 
A large MC sample of BB events is used to determine 
possible backgrounds that can contaminate our signal re- 
gion. The decay B° -^ D+n^, D+ ^ K~tt+tt+ includ- 
ing D+ -^ K*°n+ and K*°{U30)'k+, where a pion is 
misidentified as a kaon, poses a particular challenge for 
the B^ — ?■ DfiT^ channel. This decay mode has a large 
branching fraction; its reconstructed invariant mass spec- 
trum peaks near the Df peak while its AE distribution 
is shifted by about 70 MeV from zero. The B° -^ D+tt^ 
background is more prominent in Df — > K*^K^ com- 
pared to Df — > (/)7r+ because of the wider invariant- 
mass requirement. To suppress this background, we re- 
ject event candidates that are consistent with the D^ — >■ 
K^TT^TT^ mass hypothesis within 16 MeV/c^ (^ 3ct) when 
the two same-sign particles in the D'^ candidate are as- 
signed to be pions. For the J3+ — >■ K';^K^ mode there is 
a similar background from B'^ — > D^tt^ , Z?+ -^ K^Itt^. 
Here we exclude candidates consistent within 20 MeV/c^ 
with the Z?+ — > K'^'^'^ mass hypothesis. The chan- 
nel i3° — > DtK^ has a similar reflection background 



from B" -^ D+K-, D^ 



-^ R- 



T+. We apply the 
same rejection criteria, as in B^ — > Dfn^, to the invari- 
ant mass of the X^tt+tt^ system. Our invariant-mass 
veto requirements also reduce a similar background from 
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D+TT° 



Another BB background arises from charmless decays 
such as B° -^ K'^K-'K+, K^K+R-, K*°K+K-, and 
(j)K~TT^. These events peak at AE = as the final 
state is the same as signal, but have a broad nonpeak- 
ing Df mass distribution due to the absence of a Df 
in the final state. Finally, there is a crossfeed contribu- 
tion from 5° -^ DfK^ [B° -^ Z^+tt") to S" -^ D+tt" 
(_B" — >■ DfK~) due to a kaon (pion) faking a pion (kaon), 
which also needs to be considered. 

To determine the branching fractions of B'-* -^ Dfir^ 
and B'^ — )■ DfK~, we perform an unbinned extended 
maximum likelihood fit to the candidate events found in 
the selected regions of Afbc, AE, and M^+ (described 
above). The probability density functions (PDFs) are 
functions of AE and Af^+ . The extended likelihood func- 
tion is 



E^.' 



Nl 



N 



(2) 



where Yjm is the yield of event category j for Df decay 
mode m ((^7r+, K*^K^, or K'^K^), N is the total number 
of candidate events in three Df modes, and Vjm{xi) is 
the PDF evaluated for the variables x = {AE,Mjj+) 
measured for event i. To constrain the three D'^ modes 
to have a common branching fraction, we express the 
signal {j = 1) yield as 



Yir 



NsbBB„ 



(3) 



i=l j 



where Ng-^ is the number of BB events, B is the branch- 
ing fraction of B" -^ D+tt- (or, B° -> D+K-\B^ is 
the branching fraction of the D+ decay mode m [l5[ , and 
e„i is the detection efficiency of the corresponding decay 
mode. Finally to account for crossfeeds between the two 
signal channels, they are simultaneously fitted, with the 
i3° -> D+R- signal yield in the correctly reconstructed 
sample determining the normalization of the crossfeed in 
the B'^ — 7> Dfir^ fit region, and vice versa. 

There are four PDF components, each denoting an 
event category, for the Df decay modes considered: 
signal, crossfeed, combinatorial, and charmless back- 
grounds. The signal AE PDF shape is modeled with 
the sum of two Gaussian functions with a common ra- 
tio of the narrow component to the total for the three 
Df modes. We parametrize the signal M^+ distribution 
using the sum of two Gaussians with a common mean 
and ratio of areas for all the Df modes, and use the 
same PDF for both the DfTr~ and DfK~ channels. We 
use an asymmetric Gaussian to model the AE distribu- 
tion of the DfiT- {DfR-} crossfeed, that contributes to 
the B° -^ D+R- (B° -^ D+TT-) signal. Combinatorial 
background arises when a random track is combined with 
a correctly reconstructed or misreconstructed Df candi- 
date. This background is mostly from generic BB and 
continuum qq processes. To model misreconstructed Df 
candidates we use a linear function to describe the Af^+ 
distribution, while the signal Af^+ PDF shape is used for 
combinatorial background that contains correctly recon- 
structed Df candidates. The AE distribution in both 
cases is parametrized with a linear function. Charmless 
background events are characterized by a linear M^+ dis- 
tribution and a peaking AE, which is modeled with the 
signal PDF shape. For both signal and background PDF 
parametrizations, we obtain shape parameters from the 
corresponding MC samples. 

We calibrate various PDF shape parameters obtained 
from MC events using a large-statistics control sample 
of fiO -^ D+TT-; D+ ^ <^7r+, R*°R+, and R°R^- The 
peak positions and widths are adjusted based on the dif- 
ference between data and MC simulations observed in the 
control channel. We find the measured branching frac- 
tion of the control sample is in agreement with the cur- 
rent world-average value. We also cross-check our anal- 
ysis procedure by applying it to a data sample enriched 
with the Cabibbo suppressed decay 5° — > D+R~ , where 



the D+ decays to (/)7r+, K*^K+, and K^K^ . The mea- 
sured branching fraction is found to be consistent with 
the world-average value. 

Our fit in total has 32 free parameters. They are the 
branching fractions of both the decay channels (2); the 
yields of the charmless background (6), the combinatorial 
background with correctly reconstructed D'^ candidates 
(6), and the pure combinatorial background (6); and the 
slopes of the linear functions representing the nonpeak- 
ing Ai? shape (6), and the nonpeaking Mj^+ shape (6). 
Figure [5] shows results of the simultaneous fit for both 
the signal channels, projected onto A£^ and M^+. For 
/S.E {Mjj+) projections we apply the M^+ {AE) signal 
region requirement, as described earlier. In Table H] we 
summarize the fit results. The signal significance is calcu- 
lated as y^—2 ln(£o/'C,nax), where £max is the maximum 
likelihood value for the nominal data fit, and £o is the 
corresponding value with the signal yield fixed to zero. 
Including systematic errors (described below), which im- 
pact only the signal yield, into the statistical likelihood 
curve, through a Gaussian convolution, we determine the 
significance to be 8.0 and 9.2 standard deviations for 
B° -^ Df-K' and B" -^ D+K', respectively. 

Systematic uncertainties that affect our measurement 
are summarized in Table |TT1 The dominant one is the er- 
ror on the current world-average values of the Df decay 
branching fractions [15|. The remaining sources of sys- 
tematic error are the fixed PDF shapes, for which we vary 
the correction factors (applied to the peak positions and 
widths) in accordance with their errors obtained from the 
control sample i?" — )■ _D+7r~; MC statistics; the efficien- 
cies of tracking, FID, and K'^ reconstruction; the error 
on Ngg, assuming equal production of B'^B'^ and B^B^ 
pairs at the T(4S'); requirements on TZ, evaluated using 
the control sample; and the fit bias. We estimate the sys- 
tematic error due to fit bias as a linear sum of the bias 
itself and the statistical error on it, using ensembles of 
simulated experiments. 

We obtain the branching fractions B [B^ — > D+tt^) = 
(1.99 ±0.26 ±0.18) X 10-5 and B{W^D+K-) = 
(1.91 ± 0.24 ± 0.17) X 10"^, where the uncertainties arc 
statistical and systematic, respectively. These results are 
consistent with, and constitute a significant improvement 
over, our previous results [lO|. Using our measurement 
of B^ — )• DfTT~ in conjunction with the value of Cabibbo 
angle ^, tanOc = 0.2314 ±0.0021, the lattice QCD cal- 
culation of /dJId = 1.164± 0.011 [i3|, and the branch- 
ing fraction ^3(5° ^ D-it+) = (2.68 ± 0.13) x lO'^ [H], 
we obtain Rd^ = [1.71 ± O.ll(stat) ± 0.09(syst) ± 
0.02(theo)]%, where the last term accounts for the theory 
uncertainty in the fo/fD^ estimation. Uncertainties due 
to other possible SU(3) breaking effects '20|, which are 
of order (10-15)%, are not included in the quoted theory 
error. This constitutes the most precise measurement of 
Ro-n to date. The measured value of B (S" -> DfK~) 
can be understood in terms of a pure Il^-exchange contri- 



bution, which is in agreement with our recent measure- 
ment of B" -)■ D*+K- ^|. 

To conclude, using a data sample of 657 x 10^ BB 
pairs collected by Belle, we report the most precise mea- 
surement of branching fractions for the i?" — >■ DfTT~ and 
i?° — )■ DfK~ decays. This improves the precision of the 
parameter Rdtt, and thus will also improve determina- 
tions of the UT angle 03 from CP violation measurements 
in B'^ — >■ D^TT^ decays. One can use the S" — ?► Dfn^ 
result to calculate the CKM matrix element \Vub\ fol- 
lowing the prescription laid out in Ref. [7|. Our results 
supersede the previous Belle measurement |10| . 
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FIG. 2: (color online). Projections of the simultaneous fit to (a,c) B" — >■ D^tt" and (b,d) B" — >■ DfK~. (a,b) correspond to 
AE and (c,d) are the D^ mass distributions. Points with error bars show the data, the blue solid curves are the total fit result, 
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curves are the combinatorial background, and the blue dashed curves correspond to the charmless B background. 
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